Approximation of analytic functions in Korobov spaces 
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Abstract 



We study multivariate ^-approximation for a weighted Korobov space of ana- 
lytic periodic functions for which the Fourier coefficients decay exponentially fast. 
The weights are defined, in particular, in terms of two sequences a = {a,-} and 
b = {bj} of numbers no less than one. Let e i2_app ' A (ra, s) be the minimal worst- 
case error of all algorithms that use n information functionals from the class A in 
the s-variate case. We consider two classes A: the class A al1 consists of all linear 
functionals and the class A std consists of only function evaluations. 
43 " We study (EXP) exponential convergence. This means that 

e L2 - app ' A (n, s) < C{s) q WCiM) p(s) for all n, s G N 

where q G (0, 1), and C, C\,p : N — > (0,oo). If we can take p(s) = p > for all s 
" then we speak of (UEXP) uniform exponential convergence. We also study EXP and 

UEXP with (WT) weak, (PT) polynomial and (SPT) strong polynomial tractability. 
These concepts are defined as follows. Let n(e, s) be the minimal n for which 
e i2 " app ' A (n,s) < e. Then WT holds iff lim s+log £ -i_ s>00 (log n(e, s))/(s+log e" 1 ) = 0, 
PT holds iff there are c, t\,T2 such that n(e,s) < cs ri (l + log e _1 ) T2 for all s and 
e G (0, 1), and finally SPT holds iff the last estimate holds for t\ = 0. The infimum 
of T2 for which SPT holds is called the exponent r* of SPT. We prove that the 
results are the same for both classes A, and: 



• EXP holds for any a, b and uj. 

• UEXP holds iff B := Y%Li 1 /bj < oo and the largest p is 1/B. 

• WT+EXP holds iff lim^ a, = oo. 

• WT+UEXP holds iff B < oo and lira,- aj = oo. 

• The notions of PT and SPT with EXP or UEXP are equivalent, and hold iff 
B < oo and a* := lim inf j->oo (l°g a j ) 1 3 > 0. Then 

max(5,(log3)/a*) < r* < B + (log 3)/a*, 

and r* = B for a* = oo. 
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1 Introduction 



We study approximation of s-variate functions defined on the unit cube [0, l] s with the 
worst-case error measured in the L 2 norm. Multivariate approximation is a problem 
that has been studied in a vast number of papers from many different perspectives. We 
consider analytic periodic functions belonging to a weighted Korobov space. We present 
necessary and sufficient conditions on the decay of the Fourier coefficients under which 
we can achieve exponential and uniform exponential convergence with various notions of 
tractability. 

We approximate functions by algorithms that use n information evaluations. We either 
allow information evaluations from the class A al1 of all continuous linear functionals or 
from the class A std of standard information which consists of only function evaluations. 

For large s, it is important to study how the errors of algorithms depend not only 
on n but also on s. The information complexity n L2 - app ' A (s, s) is the minimal number n 
for which there exists an algorithm using n information evaluations from the class A G 
{A al1 , A std } with an error at most e in the s-variate case. The information complexity is 
proportional to the minimal cost of computing an e- approximation since linear algorithms 
are optimal and their cost is proportional to ri L2_app ' A (£, s). 

We would like to control how ri L2 ~ app ' A (£, s) depends on e~ l and s. In the standard 
study of tractability, see [TJEIIH], weak tractability means that n L2_app ' A (e, s) is not ex- 
ponentially dependent on e~ l and s. Furthermore, polynomial tractability means that 
n i2 - app ' A (£, s) is polynomially bounded by C s q e p for some C, q and p independent of 
£ G (0, 1) and s G N. If q = then we have strong polynomial tractability. 

Typically, n L ' 2 ^ app ' A (e, s) is polynomially dependent on e~ l and s for weighted classes 
of smooth functions. The notion of weighted function classes means that the successive 
variables and groups of variables are moderated by certain weights. For sufficiently fast 
decaying weights, the information complexity depends at most polynomially on s, and we 
obtain polynomial tractability, or even strong polynomial tractability. 

These notions of tractability are suitable for problems for which smoothness of func- 
tions is finite. This means that functions are different iable only finitely many times. 
Then the minimal errors of algorithms enjoy polynomial convergence and are bounded 
by C(s) n~ T , for some positive C(s) which depends only on s and some positive r which 
depends on the smoothness of functions. For many classes of such functions we know 
the largest r which grows with increasing smoothness and decreasing weights. Further- 
more, weak tractability holds if log C(s) = o(s), whereas polynomial tractability holds if 
C(s) is polynomially dependent on s, and strong polynomial tractability holds if C(s) is 
uniformly bounded in s. 

It seems to us that the case of analytic or infinitely many times different iable functions 
is also of interest. For such classes of functions we would like to replace polynomial 
convergence by exponential convergence, and study the same notions of tractability in 
terms of (1 + log£ _1 ,s) instead of More precisely, let e L2 ~ app,A (n, s) be the 

minimal worst-case error among all algorithms that use n information evaluations from a 
permissible class A in the s-variate case. By exponential convergence of the nth minimal 
approximation error we mean that 

e L ' 2 - app ' A (n,s) <C(s)q {n/Cl{s))P(s) for all n, s G N. 
Here, q G (0, 1) is independent of s, whereas C, Ci, and p are allowed to be dependent on s. 
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We speak of uniform exponential convergence if p can be replaced by a positive number 
independent of s. A priori it is not obvious what we should require about C(s),Ci(s) and 
p(s) although, clearly, the smaller C(s) and Ci(s) the better, and we would like to have 
p(s) as large as possible. Obviously, if we do not care about the dependence on s then 
the mere existence of C(s), Ci(s) and p(s) is enough. 
The last bound on e L2_app,A (n, s) yields 

n L2 - app > A (e,s) < 

Exponential convergence implies that asymptotically with respect to e tending to zero, we 
need C^log 1 ^- 1 e~ l ) information evaluations to compute an e-approximation to functions 
from the Korobov space. (Throughout the paper log means the natural logarithm and 
log r x means [logx] r .) 

Tractability with exponential or uniform exponential convergence means that we would 
like to replace e" 1 by 1 + log e" 1 and guarantee the same properties on n i2_app ' A (£:, s) as 
for the standard case. This means that (WT) weak tractability holds iff 

log n L2 - app > A (e,s) „ 
hm = 0, 

s+log e-Uoo S + log £ 1 

whereas (PT) polynomial tractability holds iff there are non-negative numbers c,t±,T2 
such that 

n L2 ~ app > A (e, s) < cs Tl (l + log e" 1 ) 7 " 2 for all s6N,e6 (0, 1). 

If ti = in the last bound we speak of (SPT) strong polynomial tractability, and then r* 
being the infimum of T2 is called the exponent of SPT. 

For instance, uniform exponential convergence implies weak tractability if 

C(s) = exp (exp (o(s))) and C±(s) = exp(o(s)) as s — > 00. 

These conditions are rather weak since C(s) can be almost doubly exponential and C\(s) 
almost exponential in s. 

Furthermore, uniform exponential convergence implies polynomial tractability if for 
some non-negative T]i and 772 we have 

C{s) = exp {0{s m )) and d(s) = 0{s m ) as s -)> 00. 

If rji = 772 = then we have strong polynomial tractability. 

Uniform exponential convergence with weak, polynomial and strong polynomial tract- 
ability was studied in the papers [2] and [1] for multivariate integration in weighted 
Korobov spaces with exponentially fast decaying Fourier coefficients. However, the notion 
of weak tractability was defined differently in a more demanding way, see Section M for 
more details. In the current paper, we deal with multivariate approximation in the worst- 
case setting for the same class of functions. We study exponential and uniform exponential 
convergence and various notions of tractability defined as above. 

We find it interesting that all results presented in this paper are exactly the same for 
both classes A al1 and A std . This is surprising since the class A std is much smaller than 
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Cx(s) 



log C(s) + loge" 
logg -1 



1 \ i/p(») 



for all s£N and e e (0, 1) 



the class A al1 . This is very good news since usually in the computational practice we can 
only use function values, i.e., the class A std . Furthermore, all our results are constructive 
or semi-constructivd^. That is, we provide algorithms that use only function values and 
for which we achieve exponential and uniform exponential convergence with WT, PT or 
SPT. The sample points used by these algorithms are from regular grids with varying 
mesh-sizes for successive variables. Such grids were also successfully used for multivariate 
integration in the previous papers [2] and [I]. 

For the Korobov class of functions / considered here, the decay of the Fourier coef- 
ficients f(h) is defined by two sequences a = {aj} and b = {bj}, and by a parameter 
oj G (0,1). Here a and b are two sequences of real numbers bounded below by 1, see 
Section [2] for further details. We assume that 



1 < oo, 



where 

u h = w^=i a ' |fcj|kj for all h = (h u h 2 ,. . . ,h 8 ) G Z s . 

We study for which (a, b, oj) we have exponential and uniform exponential convergence 
without or with various notions of tractability. It turns out that u only effects the factors 
in our estimates. These factors go to infinity as u tends to one. 

We are going to show that exponential convergence holds for any choice of a and b, 
whereas uniform exponential convergence holds iff 



oo 1 



< OO, 



3=1 1 



independently of a. Furthermore, the largest rate p(s) for exponential convergence is 
1/B(s), where 



B(s) 



y- 

j=l J 



and for uniform exponential convergence the largest rate p is \ jB. 

We prove that (WT+EXP) weak tractability with exponential convergence holds iff 



lim a 3 - = oo, 

j -too 



and (WT+UEXP) weak tractability with uniform exponential convergence holds iff 

B < oo and lim aj = oo. 

J-K50 

The notions of polynomial and strong polynomial tractability with exponential or uniform 
exponential convergence are equivalent. Furthermore, the strongest notion of tractability, 
namely strong polynomial tractability with uniform exponential convergence, holds iff 

lOg Oj ' 

B < oo and a* = liminf — — > 0, 

j^fOO j 



1 Semi-construction is only used for the class A std when we want to achieve WT with UEXP, see 
Section [ 
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and then the exponent r* of SPT satisfies 



, D log 3\ log 3 
max B, < t < B H . 



a j a 

We comment on the assumption that a* > 0. This means that the aj are exponentially 
large in j for large j. Indeed, a* > implies that for any 5 G (0, a*) there is j$ such that 

aj > exp(5j) for all j >j* s . (1) 

Obviously, it may happen that a* = oo. Then we know the exponent of SPT exactly, 

r* = B. 

Note that this happens if, for instance, aj > exp(abj) for large j and for some a > 0. 
Indeed, then 

a > hmmf — — = oo, 

since B < oo implies that liminf^oo bj/j = oo. 

The rest of the paper is structured as follows. We give detailed information on the 
Korobov space in Section [2j and on /^-approximation and tractability in Section |3j Our 
main results are summarized in Section HI The proofs for the class A al1 are in Section [6] us- 
ing preliminary observations from Section |5j The proofs for the class A std are in Section [8] 
using preliminary observations from Section [71 In Section M we compare the approxi- 
mation problem considered in this paper with the integration problem considered in [1]. 
Analyticity of functions from the Korobov space considered in this paper is shown in 
Section [TUJ 

2 The Korobov space H(K S a 

The Korobov space H(K Sia b) discussed in this section is a Hilbert space with a repro- 
ducing kernel. For general information on reproducing kernel Hilbert spaces we refer 
to [I]. 

Let a = {a>j}j>i and b = {bj}j>\ be two sequences of real positive weights such that 
bj > 1 and aj > 1 for all j = 1, 2, . . . . (2) 
Throughout the paper we assume, without loss of generality, that 

«1 < 0-2 < < . . . . 

Fix bj e (0, 1). Denote 

u h = w £?=i a ^ hj for all h = (h u h 2 , . . . , h s ) G If. 

We consider a Korobov space of complex- valued one-periodic functions defined on [0, l] s 
with a reproducing kernel of the form 

K s ,a,b( x > V)= S ^ UJ h exp(27rift ■ (x - y)) for all x,y G [0, 1} S , 



with the usual dot product 

s 

h ■ {x - y) = ^2h j (x j - yj), 
i=i 

where h* the jth components of the vectors h, x, y, respectively, and i = y—l. 

The kernel K s a b is well defined since 

s / oo \ 

\K s , a , b (x, y) I < K., a>b (x, x) = 11 1 + 2 ^ ^ < °°- (3) 
The last series is indeed finite since 

oo oo 

/i=i /i=i 

The Korobov space with reproducing kernel K s ^ ab is a reproducing kernel Hilbert 
space and is denoted by H(K s ^ b ). We suppress the dependence on u in the notation 
since u will be fixed throughout the paper and a and b will be varied. 

Clearly, functions from H(K sab ) are infinitely many times differentiable, see They 
are also analytic as shown in Section [TTJ 

For / G H(K s ^ b ) we have 

f(x) = E f( h ) exp(2vrih ■ x) for all a; G [0, l] s , 

where f(h) = Jj Q ^ /(a;) exp(— 2nih ■ x) dx is the hth Fourier coefficient. The inner 
product of / and g from H(K s ^ b ) is given by 



(f>9) H (K Sta>b ) = E f( h )9( h ) u h 



and the norm of / from H(K sab ) by 



Define the functions 

e h (a:) = exp(27ri h ■ x) u][ 2 for all x G [0, (4) 

Then {eh}hez a is a complete orthonormal basis of the Korobov space H(K sab ). 

Integration of functions from H{K s ^ a ^) was already considered in [I] and, in the case 
dj = bj = 1 for all j G N, also in [2]. In this paper we consider the problem of multivariate 
approximation in the L 2 norm which we shortly call L 2 -approximation. 
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3 ^-approximation 

In this section we consider L 2 -approximation of functions from H(K sab ). This problem 
is defined as an approximation of the embedding from the Korobov space H(K sab ) to 
the space L 2 ([0, 1} S ), i.e., 

EMB S : H(K s ^ b ) -> L 2 ([0, 1} S ) given by EMB S (/) = /. 

Without loss of generality, see, e.g., [ID], we approximate EMB S by linear algo- 
rithms A n<s of the form 

n 

A n Af) = a k L k(f) for / G H{K s , a , b ), (5) 

k=l 

where each a k is a function from Z/ 2 ([0, and each L k is a continuous linear functional 
defined on H(K sab ) from a permissible class A of information. We consider two classes: 

• A = A al1 , the class of all continuous linear functionals defined on H(K sab ). Since 
H(K a , a ,b) is a Hilbert space then for every L k G A al1 there exists a function f k from 
H{K S!aib ) such that L fc (/) = (f,fk) H{Ksab ) for all / G H(K Stflib ). 

• A = A std , the class of standard information consisting only of function evaluations. 
That is, L k G A std iff there exists x k G [0, 1} S such that L k (f) = f(x k ) for all 

feH(K 8>atb ). 

Since H(K Sya b ) is a reproducing kernel Hilbert space, function evaluations are contin- 
uous linear functionals and therefore A std C A al1 . More precisely, 

Lk{f) = f{x k ) = (f,K Sjajb (-,x k )) H(Ksab) and \\L k \\ = \\K s ^ b \\ H{Ks ^ b) = K l J^ b {x k ,x k ). 
The worst-case error of the algorithm A n>s is defined as 

^-v? mK ) t AniS ) : = sup 11/ - A„, s (/)|| . 

Let e L2_app,A (n, s) be the nth minimal worst-case error, 

e L 2 -app,A (r2;S) = inf e ^VV {H{Ks ab))An s ^ 

where the infimum is taken over all linear algorithms A UjS using information from the 
class A. For n = 0we simply approximate / by zero, and the initial error is 

e L 2 -app,A (0; s) = ||EMBs|| = gup [[/[I = I. 

This means that L 2 -approximation is well normalized for all s G N. 

We study exponential convergence in this paper. Suppose first that s G N is fixed. 
Then we hope that everyone would agree that exponential convergence for e i2 " app ' A (n,s) 
means that there exist functions q : N — > (0, 1) and p, C : N — > (0, 00) such that 

e i2 - app ' A (n, s) < C{s) q{s) nPW for all neN. 
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Obviously, the functions q(-) and p(-) are not uniquely defined. For instance, we can 
take an arbitrary number q G (0, 1), define the function C\ as 

i \ i/p(s) 

cm ' log q 



log q(s) 
and then 

C(s)q(s) nPM =C(s)q^ c ^ P<s \ 

We prefer to work with the latter bound which was already considered in [I] for multi- 
variate integration. 

We say that we achieve exponential convergence for e L2_app,A (n, s) if there exist a 
number q G (0, 1) and functions p, C, C\ : N — > (0, oo) such that 

e L 2 -ap P ,A^ g ) < c(s) g (n/ClW)pW for all n G N. (6) 

If dBJ) holds we would like to find the largest possible rate p(s) of exponential convergence 
defined as 

p*(s) = sup{p(s) : p(s) satisfies (jHJ) }• (7) 

We say that we achieve uniform exponential convergence for e i2_app,A (n, s) if the func- 
tion p in (E]) can be taken as a constant function, i.e., p(s) = p > for all s G N. Similarly, 
let 

p* = sup{ p : p(s) = p > satisfies (jHJ) for all s G N } 
denote the largest rate of uniform exponential convergence. 

For e G (0, 1), s G N, and A G {A al1 , A std }, the information complexity is defined as 

n L2 ~ app ' A (e, s) := min {n : e L2 " app ' A (n, s) < e} . 

Hence, n i2_app ' A (£, s) is the minimal number of information evaluations from A which 
is required to reduce the initial error eQ 2 ~ app , which is one in our case, by a factor of 
s G (0,1). Clearly 

n i2 - app ' AStd (£,s) > n L2 - app < Aall (e,s). 

We are ready to define tractability concepts similarly as in [2] and [1] . We stress again 
that these concepts correspond to the standard concepts of tractability with e _1 replaced 
by 1 + log e -1 . We say that we have: 

• Weak Tractability (WT) if 

log n L2 " app ' A (e,s) 
hm = 0. 

s+log e-^oo S + log E 1 

Here we set log = by convention. 

• Polynomial Tractability (PT) if there exist non-negative numbers c,t\,ti such that 

n i2 - app>A (£,s) < cs Tl (1 + log e^Y 2 for all s G N, e G (0,1). 
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• Strong Polynomial Tractability (SPT) if there exist non-negative numbers c and r 
such that 

n L 2 - app ,A^ s )<c(l + log e' 1 ) T for all s G N, e G (0, 1). 

The exponent r* of strong polynomial tractability is defined as the infimum of r for 
which strong polynomial tractability holds. 

A few comments of these notions are in order. As in [2] , we note that if ([6]) holds then 

n La ~ app ' A (£, s) < 

Furthermore, if ([8]) holds then 

e i2 - app ' A (n + 1, s) < C(s)q {n/Cl{s))Pis) for all s,n G N. 

This means that ([6]) and ([8]) are practically equivalent. Note that l/p(s) determines the 
power of log e~ 1 in the information complexity, whereas log q~ x effects only the multiplier 
of log 1 ^ 5 -* e -1 . From this point of view, p(s) is more important than q. That is why we 
would like to have ()6]) with the largest possible p(s). We shall see how to find such p(s) 
for the parameters (a, b, u) of the weighted Korobov space. 

Exponential convergence implies that asymptotically, with respect to e tending to 
zero, we need 0(log 1//p ^- ) information evaluations to compute an e- approximation to 
functions from the Korobov space. However, it is not clear how long we have to wait to 
see this nice asymptotic behavior especially for large s. This, of course, depends on how 
C(s),Ci(s) and p(s) depend on s. This is the subject of tractability which is extensively 
studied in many papers. So far tractability has been studied in terms of s and e _1 . The 
current state of the art on tractability can be found in [TJ O |9]. In this paper we follow 
the approach of [2J and [4J and we study tractability in terms of s and 1 + loge" 1 . In 
particular, weak tractability means that we rule out the cases for which n L2 ~ app ' A (e, s) 
depends exponentially on s and log e™ 1 . 

For instance, assume that (Q holds. Then uniform exponential convergence implies 
weak tractability if 

C(s) = exp (exp (o(s))) and Ci(s) = exp(o(s)) as s — > oo. 

These conditions are rather weak since C(s) can be almost doubly exponential and C\(s) 
almost exponential in s. The definition of polynomial (and strong polynomial) tractability 
implies that we have uniform exponential convergence with C(s) = e (where e denotes 
exp(l)), q = 1/e, Ci(s) = cs Tl andp = l/r 2 . For strong polynomial tractability C\(s) = c 
and r* < 1/p*. 

If (jHJ) holds then we have polynomial tractability if p := mf s p(s) > and there exist 
non-negative numbers A, A\ and 7], r\\ such that 

C(s) < exp (As 11 ) and d(s) < A 1 s Vl for all s G N. 

The condition on C(s) seems to be quite weak since even for singly exponential C(s) 
we have polynomial tractability. Then r± — T)\ + T)/p and r 2 = 1/p. Strong polynomial 
tractability holds if C(s) and C*i(s) are uniformly bounded in s, and then t* < 1/p. 
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CAs) 



logC(s) + log£" 

log q^ 1 



l/p(s) 



for all s6N and e G (0,1). (8) 



4 The main results 



We first present the main results of this paper. We will be using the following notational 
abbreviations 

EXP UEXP WT PT SPT 
WT+EXP PT+EXP SPT+EXP 
WT+UEXP PT+UEXP SPT+UEXP 

to denote exponential and uniform exponential convergence, weak, polynomial and strong 
polynomial tractability, as well as weak, polynomial and strong polynomial tractability 
with exponential or uniform exponential convergence. We want to find relations between 
these concepts as well as necessary and sufficient conditions on a and b for which these 
concepts hold. As we shall see, many of these concepts are equivalent. 

Theorem 1 Consider L 2 - approximation defined over the Korobov space with kernel K s a b 
with arbitrary sequences a and b satisfying (|2J). The following results hold for both classes 
A al1 and A std . 

1 EXP holds for arbitrary a and b and 

s 1 

p*{s) = l/B(s) with B(s) := Jj— . 

j=l 3 

This implies that 

WT <£> WT+EXP, PT <£> PT+EXP, SPT & SPT+EXP. 



2 UEXP holds iff a is an arbitrary sequence and b such that 

oo 1 



< OO. 

3=1 JJ 

If so then p* = 1/B and 

WT & WT+UEXP, PT <£> PT+UEXP, SPT & SPT+UEXP. 

3 Polynomial (and, of course, strong polynomial) tractability implies uniform expo- 
nential convergence, PT =>■ UEXP, i.e., 

PT <£> PT+UEXP, SPT & SPT+UEXP. 

4 We have 

WT lim Oj = oo, 

WT+UEXP & B < oo and lim a,- = oo. 

5 The following notions are equivalent: 

PT <=> PT+EXP & PT+UEXP SPT <=> SPT+EXP SPT+UEXP. 
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6 SPT+UEXP holds iffbj 's are summable and aj 's are exponentially large in j, i.e., 



oo 



log aj 
3 




— < oo and a* := liminf 



> 0. 



Then the exponent r* of SPT satisfies 



max I B, 



log 3 



a 



) 



<r* <B + 



log 3 



a* 



In particular, if a* = oo then r* = B. 

We comment on Theorem [U We already expressed our surprise in the introduction 
that the results are the same for both classes A std and A al1 , although the class A std is much 
smaller than the class A al1 . However, the proofs for both classes are different. We also 
stress that the results are constructive. The corresponding algorithms can be found in 
Section for the class A al1 and in Section [S] for the class A std . 

Point Q] tells us that we always have exponential convergence and the best rate is 
p*(s) = l/B(s). Note that p*(s) decays with s, and if B(s) goes to infinity then the rate 
decays to zero. The smallest rate is for bj = 1 for all j > 1, for which p* (s) = 1/ s. Clearly, 
all tractability notions with or without exponential convergence are trivially equivalent. 

Point [2] addresses uniform exponential convergence which holds iff bj lj s are summable, 
i.e., when B < oo. Then the best rate of uniform exponential convergence is p* = 1/B. 
Obviously, for large B this rate is poor. We stress that uniform exponential convergence 
holds independently of a. Similarly as before, as long as B < oo, tractability notions with 
or without uniform exponential convergence are trivially equivalent. 

Point [3] states that (strong) polynomial tractability implies uniform exponential con- 
vergence, i.e., B < oo. This means that the notion of polynomial tractability is stronger 
than the notion of uniform convergence. 

Point H] addresses weak tractability which holds iff a^-'s tend to infinity. We stress 
that this holds independently of b and independently of the rate of convergence of a to 
infinity. We have weak tractability with uniform convergence if additionally B < oo. 
Hence for linx, aj = oo and B = oo, weak tractability holds without uniform exponential 
convergence. 

Point [5] states that, in particular, the notions of polynomial tractability and strong 
polynomial tractability with uniform exponential convergence are equivalent. 

Point [H] presents necessary and sufficient conditions on strong polynomial tractability 
with uniform exponential convergence. We must assume that B < oo and a* > 0. The 
last condition means that a/s are exponentially large in j for large j. We only know 
bounds of the exponent r* of strong polynomial tractability. Note that for large B or 
small a* the exponent r* is large. On the other hand, r* is not large if B is not large 
and a* is not small. We stress that B can be sufficiently small if all bj are sufficiently 
large, whereas a* can be sufficiently large if aj are large enough. In fact, we may even 
have a* = oo. This holds if aj goes to infinity faster than C- 7 for any C > 1. We already 
noticed in the introduction that this holds, for example, if aj > exp(5bj) for large j and 
for some 5 > 0. For a* = oo we know the exponent of SPT exactly, 



r* = B. 
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5 Preliminaries for the class A al1 

The information complexity is known for the class A al1 , see, e.g., [XQ1 Chapter 4, Sec- 
tion 5.8]). It depends on the eigenpairs of the operator 

W s = EMB: EMB s : H{K Sta>b ) H(K Sta>b ), 

which in our case is given by 

with given by fll]). Hence, the eigenpairs of W s are {uhi e h) since 

W s e h = oo h e h = oj^U a ^ e h for all h G Z s . 

It is known that the information complexity is the number of the eigenvalues Uh of the 
operator W s which are greater than e 2 . More precisely, for a real M define the set 

A{s,M) := {heZ s : oo^ < M) 

= {heZ s : W -Sj=i°il h il ij < M| . (9) 

Then 

n i2 - app ' AaI1 (e, s) = |^(s,^ 2 )| . (10) 
Furthermore, the optimal algorithm in the class A al1 is the truncated Fourier series 

4? } Cf)(*):= E </• e *>*(JW) e * = E /Wexp(2 7 ri^. ;E ) ) 

/ie^(s,e- 2 ) /ie^l(s,e- 2 ) 

where n = |.A(s, e -2 )|, which ensures that the worst-case error satisfies 

e L2 -^(H(K s>a>b ),A^f) <e. 

For the proof of Theorem [1] and also for the further considerations in this paper we 
need a few properties of the set A(s, M) and its cardinality. Clearly, A(s, M) = for all 
M < 1. For e G (0,1), let 

. , log 6~ 2 

x = x(e) := - > 0, 

log UJ~ L 

and 

( s 

I \ I > _ r77S V^ 1 IT lh„- 



n(x, s) :- 
Then 



i=i 



n i2 - app ' A (e,s) = \A(s,e- 2 )\=n(x,s). 
We have n(x, s) = 1 for all x G (0, ai] and 

n(x,l) = 2 \(x/ai) 1/bl ] - 1, 
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n(x,s) = n(x, s — 1) + 2 n(x — a s h bs , s — 1). 

h=l 

Clearly, n(y, s) > n(x, s) > n(x, s — 1) for all y > x > and s > 2. Note that for 
x < a s , the last sum in n(x, s) is zero and n(x, s) = n(x, s — 1). For x > oi, define 

j(x) = sup{ j G N : x > % }. 

For linij a, < oo we have = oo for large x. For linij a,j = oo, we can replace the 
supremum in j(x) by the maximum, and j(x) is finite for all x. However, j(x) tends to 
infinity with x. 

If j (x) is finite then 

n(x, s) = n(x,j(x)) for all s > j(x), 
and therefore, if j(x) < oo then 

hm l ° gn{X ' S) = 0. 

s^oo S + X 

We now prove the following lemma. 
Lemma 1 

• For x > ai + a 2 + • • • + a s we have 

n{x,s) > 3 s . 



For x > ai and for arbitrary aj G [0, 1] we have 



x 



x 



min(s,j(x)) I 

n(x,s) > Jl I 2 

min(s,j(x)) / 

n(x,s) < 2 

j=i V 

where the empty product is defined to be 1. 
• For x > ai we have 



i/bj 



' fc=i+i 



1 . 



n( 2 



min(s ,j(x)) 

1 ) < //(.r. s) < Yl [ 2 



l/^ 



- 1 



Proof. To prove the first point, let A s = { h G Z s : ^ G {-1, 0, 1} }. For h e A s we 
have 

s s 

i=i i=i 
Hence 3 s = |Ag| < n(x, s), as claimed. 
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We turn to the second point. It is easier to prove the upper bound on n(x, s). From 
the recurrence relation on n(x, s) we have 



n(x,s) < n(x, s — 1)+2 
This yields 



l/b s 



1 ] n(x, s— 1) 



l/b a 



1 nix, s—1). 



n(x,s)<H[2 

3=1 

If j > j(x), i.e., x < dj, then the factor 

1/6/ 



X 



1/bj 



1 . 



X 



-1 = 2-1-1 = 1. 



Hence, we can restrict j in the last product to min(s, j(x)) and obtain the desired upper 
bound on n(x, s). 

We turn to the lower bound on n(x, s). Note that x — a s h bs > a s x for all h £ N with 

'x(l-a s )^ 1/bs 



h < 



Hence 



n(x,s) > n(a s x, s — 1) + 2n(a s x, s — 1) 



1. 



x(l — a. 



l/bs 



x(l — a, 



rt,. 



— 1 n(a.x, s — 1). 



We now apply induction on s. For s = 1 we have 

n(x, 1) = 2 \(x/ ai ) 1/bl ] - 1 > 2 [((1 - ai)x/a!) 1/61 ] - 1, 



as claimed. Then 



nfx, s) > 2 



x(l-a s )V A 



min(sj(:r)) 

= n 



.r 



Clj 



0S« a 

i 1/6, 

fc=j+l 



s-1 



1 - 



n a A 



\ lf>j 



k=j+l 



X 



[1 - (Xj) Yl a k 
k=j+l 



- 1 



as claimed. This completes the proof of the second point. 

To prove the third point, it is enough to take aj = (j — Then for j = 1, 2, 
we have 

a nu-#-i) i 
11 ak - 



[1 - a« 



k=j+i 



jUl= j+ i k s l 



as claimed. This completes the proof of Lemma [TJ 



□ 
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6 The proof of Theorem D for A al1 

We are ready to prove Theorem [1] for the class A al1 . 

6.1 The proof of Point U 

From the second and third points of Lemma Q] with a fixed s we have 

n(x,s) = Q(x B( - s ^) as x — > oo. 
Therefore there are functions c\, ci : N — > (0, oo) such that 

ci(s) \og B{s) e~ l < n L2 - app ' Aall (e,s) < c 2 (s) log B(s) s' 1 
for e tending to zero. This implies exponential convergence since 

e L2 - app ' Aa "(n, s) < q ^/c2{s)Y' B ^ with q = exp (_!). 
Hence, p*(s) > 1/B(s). On the other hand, if we have exponential convergence (jS]) then 

n i2 - app ' Aall (e, S ) = O (log 1/p(s) 

and l/p(s) > B(s), or equivalently, p(s) < 1/B(s). Hence, p*(s) = l/B(s), as claimed in 
Point [TJ The rest in this point is clear. 

6.2 The proof of Point H 

Assume now that we have uniform exponential convergence. Then e i2_app ' Aall (n, s) < 
C(s) g( n / c ' 1 ( s )) p implies for a fixed s that 

n L2 " app ' Aall (e, s) = C(log 1/p a' 1 ) as e -> 0. 

Then B(s) <l/p for all s. Therefore B < 1/p < oo and < 1/B. On the other hand, if 
B < oo then we can set = 1/B and obtain uniform exponential convergence. Hence, 
p* > 1/B, and therefore p* = 1/B, as claimed. The rest of Point [2] is clear. 

6.3 The proof of Point [3] 

PT means that 

n i2 ~ app ' Aall (e,s) < cs T1 (1 + log e" 1 ) 7 " 2 . 

This implies that 

e L 2 -app,A all ( n ) < ^-(n/c^i) 1 /^ 

Hence, UEXP holds with p = l/r 2 . This also yields the equivalence between various 
notions of tractability with or without uniform exponential convergence. 
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6.4 The proof of Point H 

We first prove that WT implies lint, dj = oo. We use the first part of Lemma [TJ For 
5 > 0, take x = (1 + $){a\ + • - • + a s ), or equivalently 

log e' 1 = - loguT 1 = - (1 + 5)(ai + a 2 H h a s ). 

Then 

log n(x, s) s log 3 log 3 



s + log e 1 s + | x log 1 + |(1 + 5)?/ s log uj- 1 ' 

where 

ai + a 2 + ha s 

Vs = • 

s 

WT implies that lim s z s = 0. This can hold only if lim s y s = oo which implies that 
lira, a,j = oo, as claimed. 

Next, we need to prove that linij aj = oo implies WT. The eigenvalues of W s are Uh 
for all li 6 Z s . Let the ordered eigenvalues of W s be A S; „ for n 6 N with A s .i > A s 2 > 
A s , 3 > — Obviously {A Si „}„ e N = {uh}h£Z a - Therefore for any rj G (0, 1) we have 



3=1 heZ 3 3=1 \ h=l 

Note that 



This proves that 



h=l h=l 



TT3 1 (l + 2w" a */(l-w'"«)) 
A s , n < (11) 



Since n L2 ~ app ' Aall (e, s) = min{n : A Sj „ + i < e 2 } we conclude that 

rn 1 (1+2^/(1 -wi**)) 

n i2 - app - A (e,s) < -. 

Using log(l + x) < x for x > 0, this yields 

s 

log n i2 - app ' Aall (e,s) < 277 log s~ l + 2 ^c i; 

i=i 

where 



1 - u;^ 

Note that linx, a, = oo implies that linx,- Cj = 0, and lim s Y^j=i c 3 -/s = 0. Therefore 

log n L2 - app ' Aa "(e,s) 
hmsup — <2?7. 

s+log e-^oo s + log £ 

Since 77 can be arbitrarily small this proves that 

log n i2 - app ' Aall (e,s) „ 
hm = 0. 

s+log e -i->oo S + log E 1 

Hence, WT holds for the class A al1 , as claimed. The rest in this point follows from the 
previous results. This completes the proof of Point HI 
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6.5 The proof of Points [5] and [6] 



For Point [5j it is enough to prove that PT implies SPT+UEXP. This will be done by 
showing that PT implies that B < oo and a* > 0. Then we show that B < oo and a* > 
imply SPT+UEXP and obtain bounds on the exponent of SPT. 

We know that PT implies UEXP and that UEXP implies that B < oo. From the 
lower bound of Lemma [T] with x — (1 + S)(ai + ■ ■ ■ + a s ) and from PT we have 



1 + 6 



3 s < n{x,s) <Cs Tl 1 + ——(log w _1 )(ai + • • • + a s ) 



for all seN. 



Since ai < a 2 < a 3 < . . ., this yields 



s a s > ai + ■ ■ ■ + a s > 



1 + 5) log u- 1 



Cs^ 



l/r 2 



for all sGN. 



Hence, 



a 



liminf — - — - > — - — > 0, 



s— i-oo S T"2 

as needed. This also shows that r 2 > (log 3) /a*. Since this holds for all r 2 for which 
we have SPT, we conclude that the exponent r* of SPT also satisfies r* > (log3)/a*. 
Clearly, r* cannot be smaller than the reciprocal of the exponent p* of UEXP. Hence, 
t* > B. This completes this part of the proof as well as the proof of lower bounds on the 
exponent of SPT. 

Assume now that B < oo and a* e (0, oo]. From ([!]) with 5 G (0, a*) we have 



aj > exp(5j) for all j > f s . 



Then 



j(x) < max jg 



.* logo; 
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For x > ax, the upper bound on n(x, s) from Lemma [T] yields 



n(x, s) < 



< 




gmin(s ,j(x)) 



< x max 



< 3 j s x B+{log 3)/<5 . 



(12) 



Hence, SPT+UEXP holds, as claimed. Furthermore, since 6 can be arbitrarily close to 
a*, we conclude that the exponent of SPT satisfies 



r* < B + 



log 3 



where for a* = oo we have = 0. This completes the proof of Point [5] and of Point [6j 
The proof of the whole theorem for the class A al1 is now completed. □ 
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7 Preliminaries for the class A std 



We state some preliminary observations which will be needed to prove Theorem [T] for the 
class A std . Based on the definition of the set A(s, M) in (Q for M > 1, we will study 
approximating / G H(K s a b) by algorithms of the form 

A n ,s,M(f)( x ) = E {-^2f( x k)exp(-2nih-x k )) exp(2nih ■ x), (13) 
heA(s,M) V fc=l / 

where a? G [0, l] s . Note that A n>S} M is a linear algorithm as in §5§ with 

— exp(27rih • {x — £Cfc)) 



he.A(a,Af) 



and with Lf-(f) = f(xk) for deterministically chosen sample points x k G [0, l) s for 1 < 
k < n. Hence, L k G A std . The choice of M and x k will be given later. 

We first study upper bounds on the worst-case error of A n >StM . The following analysis 
is similar to that in [5]. We have 

(f - A n>S)M (f))(x) = E f(h)exp(2nih-x) 

hgA(s,M) 

+ E lf{h) - -^2f{x k )exp(-2nih ■ x k )) exp(27iih ■ x). 

heA(s,M) V k=l J 

Using Parseval's identity we obtain 

|2 



11/ — ^n,s,Af(/)||L 2 ([0,l] s ) 

= E i/>>i 2 + E 

h<£A{s,M) heA{s,M) 

= E i/>)i 2 



i n 

f{h) 2j /(**) exp(-2vri/i ■ a; fc ; 



fe=l 



+ E 

heA(s,M) 



f 1 n 

/ /(a?) exp(— 27rih ■ x)dx \. f( x k) exp(— 27ri/i • x k 



We have 



.(14) 



(15) 



For the second term in (TT4"1) . we make a specific choice for the points X\, . . . ,x n used in 
the algorithm A n ^ S) M- Namely, we take x^s from a regular grid with different mesh-sizes 
for successive variables. Such regular grids have already been studied in [21 H]. We now 
recall their definition. For s G N, a regular grid with mesh-sizes m%, . . . , m s G N is defined 
as the point set 



{(fci/roi, . . . , k s /m s ) : kj — 0, 1, . . . , rrij — 1 for all j = 1, 2, . . . , s} 
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where n = Ylj=i m j * s the cardinality of Q n<s . By Q^ a we denote the dual of G n ,si i.e., 
Q^ s = {h G If : hj = (mod rrij) for all j = 1,2,..., s}. 
We will make use of the following result whose easy proof is omitted. 
Lemma 2 Let Q nyS = {xi, . . . , x n } be defined as above. For any f G H(K s a fj) we have 



r i n 

/ f(x)dx--Y,f(x k ) 



E Tw 

he9i s \{0} 



For h G Z s define fh{x) : = /(a?) exp(— l-nxh • cc). Note that with / also fh belongs to 
H{K s , a ,b) and that fh{k) = f(h + k). From Lemma [2] we obtain 



/ fh(x)dx — y2fh(x k ) 



E AO 

lef&.\{0} 



E /e+'O 

«6^„\{0} 



< ( E + 



hli 1 1 E 



< 



*eei.\{0} 



Therefore, and using f JT4|) and ffT5]) for any / G H(K s ^b) with ||/||^(a: s a b ) < 1, we obtain 



\\f ~ A n,s,M(f)\\L 2 ([0,lY) - M 



h+ E E 



(16) 



It is easy to see that 



\£\ b < 2 b (\h + £\ b + \h\ 
for any h, i G Z and any & 6 N. For /i G A(s, M) this implies 

Using ([TBI) . ffTTl) and Lemma CD with x = (logM)/(logo; _1 ), we obtain for any / G 
with \\f\\ H{Ksab) < 1, 

||/-A niSi M(/)||i 2([0|in < ^ + M\A(s,M)\ u^^W 



^ M +M (n( i + 2 



logM 
a,- logo; -1 



1/6,- 



where 



F n := ^ 
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This means that 



[e L ^(H(K s>a>b ), A n , sM )f < — + M ( J] ( 1 + 2 



\j=l 



logM 
cij logo; -1 



(18) 



Furthermore, 



rr I , log M VA .A / /log AO 1 "' 



cijloguj 1 J J \ ^logu; -1 



< 2 s Y[ (l + log" 1 ^ or 1 ) Y[ (l + log 1 ^ ikf) . 



3=1 3=1 

Since M is assumed to be at least 1, we can bound 1 + log 1/6j M < 2M 1 /^ , and obtain 



n 1+2 



logM 
dj log w^ 1 



i/bj 



s 

< 4 s M B(s) Y[ (l + log 



where, as in the previous sections, B(s) := ^2j =1 bj ■ Plugging this into (ITHjl . we obtain 



[e i2 - app (^(^, ai6 ),A niM// )] 2 < i + M B M +1 D( S ,u;,&)F n , 



(19) 



where 



D( S ,a;,b):=4 s J](l + log^ 



8 The proof of Theorem U for A std 

We now present the proofs for the successive points of Theorem [T] for the class A std . 

8.1 The proof of Point U 

The following proposition will be helpful. 
Proposition 1 For s G N and £ G (0, 1) define 



m = max 

j=l,2,...,a 



log (l + 



2s 



log(l+r,2) 



B(s) 



log U 



-1 



where 



B(s)+2 
2 



Let Q* s be a regular grid with mesh-sizes mi,m 2 , . . . , m s given by 

s 

rrij := \jn l ^ B ^' b ^\ for j — 1, 2, . . . , s and n = Y\ m j 

3=1 
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Then for M = 2/e 2 we have 

e i2 - app (#(^,a,6), A njSjM ) <e, and n = (\og B{s) (l + e 
with the factor in the O notation independent of e^ 1 but dependent on s. 
Proof. We can write 



s / oo 
-b; \o \h; -i— r / v — -\ n-i> 



Fn= ^ u t.u^^p =-i+n 1+2 



dj 2 i (nij h) 



l&i.\{0} J=l \ h=l 

Since [xj > ^/2 for all x > 1, we have 

•Af' > {\hj\/2p m 1/B{s] for all j = 1, 2, . . . , s. 



Hence, 



^<-i+n 1+2 S 

i=i V h=i 



oo 



Since 6 3 - > 1 we further estimate 



h=l h=l L-U > 



From the definition of m we have 



^Q g (lW) fora „ , U2 , 



1 _ W m 1 / B (=)a j r l, 3 - 2s 

This proves 

< -1 + ( 1 + lQg(1 + r?2) l < -1 + exp(log(l + r/ 2 )) = 7? : 



Now, plugging this into (119]) . we obtain 



[e L2 - app (^(^ a , 5 ),A„, s , M )] 2 < -L + M B W +1 J D( S ,a;,% 2 . 



Note that 

1 



D(s, u, 6)^+2^5^+2 
Hence we are allowed to choose 

1 



M 



1 2 1 



D(s,u, b) s M+ 2 ?7 s ( s )+ 2 
which yields, inserting into (12"T|) . 



1 2 



[e L2 - app (^(fr s , a , 6 ), A n , SiM )] 2 < 2,0(5,0;, b)OTr/OT = £ 2 
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as claimed. 

It remains to verify that n is of the order stated in the proposition. Note that 

s s 
3=1 3=1 

However, as pointed out in [I], 

m = o(\og B{s) (l + tj- 1 )), 
as i] tends to zero. From this, it is easy to see that we indeed have 

m = O (\og B{s) (1 + e- 1 )), 
which concludes the proof of Proposition [TJ □ 

To show Point Q] for the class A std , we conclude from Proposition [1] that 



n 



This implies that we indeed have exponential convergence for A std for all a and b, with 
p(s) = 1/B(s), and thus p*(s) > 1/B(s). On the other hand, note that obviously 
e L2_app ' Astd (n, s) > e i2_app,Aall (n, s), hence the rate of exponential convergence for A std 
cannot be larger than for A al1 which is 1/B(s). Thus, also for the class A std we have 
p*(s) = 1/B(s). The rest of Point [T] is clear. 

8.2 The proof of Point H 

We turn to Point [2] for the class A std . Suppose first that a is an arbitrary sequence and 
that b is such that 

< OO. 

3=1 J 

Then we can replace B(s) by B in Proposition (TJ and we obtain 



n L 2 -app,A std ( £) S) = Q ( log B (! + e -l^ ; 

hence uniform exponential convergence with p* > 1/B holds. On the other hand, if we 
have uniform exponential convergence for A std , this implies uniform exponential conver- 
gence for A al1 , which in turn implies that B < oo and that p* < 1/B. The rest of Point [2] 
follows immediately. 

8.3 The proof of Point [3] 

The proof of Point [3] follows by the same arguments as for A al1 . 
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8.4 The proof of Point H 

We now prove the first part of Point H]for the class A std . Assume that WT holds for the 
class A std . Then WT also holds for the class A al1 and this implies that lim,, aj = oo, as 
claimed. 

Assume now that linx, aj = oo. We use [9j Theorem 26.18]. In particular, this theorem 
states that if the ordered eigenvalues A Sjn 's of W s satisfy 



M" 

K,n < for all nGN, (22) 



for some positive M SjT and r > ~ then there is a semi-constructive algorithmic such that 

e L ^ AStd (n + 2,s)<^^ for all neN (23) 

where C(r) is given explicitly in [HI Theorem 26.18]. However, the form of C(r) is not 
important for our consideration. 

For 7] G (0, 1), let r = 1/ (2rj) > ~. We stress that r can be arbitrarily large if we take 
sufficiently small rj. We already showed in the proof for the class A al1 , see (fTTj) . that we 
can take 



M s,r = n 0- + 2c ?) r < 00 with C i = 



Furthermore, we know that lim,- aj = oo implies that lim s 2~^j=i Cj/s = 0. 
From ([23]) we obtain 



L 2 -app,A 3td 



(e,s) < 3 + (M s , r C(r)) (1+1/(2r))/r e ^ 1+1 ^ T . 



This yields that 

log n L2 ' app ' Astd (e,s) ^ \\ \ x ,._ log M s 

s+log e _1 — >oo 



lim sup : < [ 1 H — 1 + lim sup 

s + log e^ 1 \ 2t J t v 



Since (logM SjT )/s < 2r $^ =1 c i/ s tends to zero as s — > oo, we have 

log n L2 - ap P' Astd (e,s) / 1 \ 1 

limsup — !-!_>< (i + _) _. 

s+ iog e -i^oo s + log£ 1 V 2ry r 

Since r can be arbitrarily large this proves that 

lim ° g f^- = 0. 

s+log e-i-foo S + log £ _1 

This means that WT holds for the class A std , as claimed. 

We turn to the second part of Point H] for the class A std . This point easily follows from 
the already proved facts that WT holds iff linr, aj = oo and UEXP holds iff B < oo. 



2 By semi-constructive we mean that this algorithm can be constructed after a few random selections 
of sample points, more can be found in [5]. 
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8.5 The proof of Point 

Suppose that PT holds for the class A std . Then PT holds for the class A al1 . By Point [5] for 
the class A al1 , which has already been proved, this implies SPT+UEXP for the class A al1 
which in turn implies that B < oo and a* > by Point [H]for the class A al1 . This implies 
SPT+UEXP for the class A std as will be shown in the subsequent Section 18.61 The rest 
of this point is clear. 



8.6 The proof of Point M 

The necessity of the conditions for SPT+UEXP on b and a stated in Point for the 
class A std follows from the same conditions for the class A al1 and the fact that the infor- 
mation complexity for A std cannot be smaller than for A al1 . 

To prove the sufficiency of the conditions for SPT+UEXP on b and a stated in Point [6] 
we analyze the algorithm A us m given by f[T51) . where the sample points are from the 
regular grid Q n)S with mesh-sizes 



mj = 2 



log M 

aj log uj- 1 



for all j — 1, 2, . . . , s. 



Here M > 1 and j3 G (0, 1). Note that rrij > 1 and is always an odd number. Furthermore 
rrij = 1 if dj > ((logM)/(logw- 1 )) 1//3 . Assume that a* G (0,oo]. Since for all 5 G (0, a*) 
we have 



see ((I]), we conclude that 



a,j > exp(5j) for all j > j* s , 



log(((lo g M)/(loga;- 1 )) 1 /0 ^ 
! " 8 J 



implies rrij = 1. 



J > 3 p,8 ■= m ax ( j s , 
From ffTBT) we have 

[e L >-^(H(K Sta>b ),A niS>M )} 2 <±- + Y Y w w- 



heA(s,M) leQlMo} 



We now estimate 



E Uh + l = 
«e^.\{0} 0^uc{i,..., s } jeu \* jG z\{o} 



E II E oo a ^ +m ^ lJ Y[ 



where we separated the cases for £j G Z\ {0} and £j = 0. We estimate the second product 
by one so that 



Y <*™ ± E II E 



*e£+ s \{o} 



0^uC{l,...,+ j'Gu W 6 Z\{0} 
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We now show that for h e A(s, M) we have \hj\ < (rrij + l)/2 for all j = 1,2, ... ,s. 
Indeed, the vector h satisfies rij=i w~~ a ^ h ^ 3 < M, and since each factor is at least one we 
have u^ a ^ h ^ 3 < M for all j, which is equivalent to 



as claimed. 

In particular, if rrij = 1 then hj = and 

y u; a ^ +m ^ bj = 2 y u a > £bj < 2 y ^ = 2ujaj < (24) 

££Z\{0} l=\ l=\ 

Let rrij > 3. Then \hj\ < (rrij + l)/2. Since both \hj\ and (rrij + l)/2 are positive 
integers, we conclude that \hj\ < (rrij + l)/2 — 1 = (rrij — l)/2 and I 7^ implies 

m ■ -I- 1 

\hj + m^\ > m 3 \i\ -\hj\> 3 

Therefore 

00 o a ■ \(m +1) /2l 6 J 

^hj+m^ ^ 2 ^ UJ a^m J + l)/2] b n b J < ^ J J (25) 

£ez\{o} 1=1 ^ U 

The inequalities (j2"lj) and fT2"5"j) can be combined as 

o, ,a,[(m, + l)/2] 6 J 



l-uo 

e&z\{o} 



Note that 



e n e ^ +mj " £ n^=- i +ri( i +^)- 

0^uc{i,..., s } jeu yez\{0} / uc{i,..., s } jeu j=i 

Consequently, 

2 bJ a J [{ m] +l)/2] b J 



4 >s <^+i^,M ) ^-i + n^+^ r 



Using log(l + x) < x we obtain 



log 



2 u a j [(m j +l)/2p 



J=l \ / J 3=1 



s 



^[(mj+l)/2]°3 

3=1 



From the definition of rrij we have aj[(rrij + l)/2] bj > a 1 - 13 (log M)/log uj 1 . Therefore 



1 V 

M 
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Since aj > 1 for j < jt s — 1 and Oj > exp(5j) for j > s we obtain 



~ . = _Vm%[K+i)/2)]' j < 1 3 1A 1 + V - ' - M 

i =1 \ 3=J},6 V 7 



M 



where 



a 



I 00 / 1 \ ex p(( 1 ~' 3 )' 5 -?)~ 1 

" ^- 1+ E(m) ' x 



3—3/3,5 



Note that for M > (7^$ we have 7 < 1. 

Using convexity we easily check that — 1 + exp(7) < (e — 1)7 for all 7 G [0,1]. Thus 
for M > we obtain 



s ( n, ,a j [(m j +l)/2] b i\ ( 9 s 



< 



-1 + exp(7) < (e — 1)7 
(e ~ 1) 
Af 



We now turn to \ A(s, M)\ which was already estimated in the proof for the class A al1 , see 
(H2D- We have 

B+(log3)/5 



\A(s,M)\ < $>(>.* 1 + 



log M 

log UJ~ X 



Therefore 



e 2 < — 

n,s - M 



1+Cw(e - 1) 3^( 1 + ^ r ) 



B+(log3)/5" 



< 



D 



0,6 



M 



where 



D PtS := sup 

x>Cp >s 



+ 



log x 
log 



B+(log3)/5 N 



< OO. 



Hence for 
we have 



M = max(C^, J D^e- 4 ) 



e„, fl = e L2 - app (i/(^ ia , f) ),A niS ,M) < e. 
We estimate the number n of function values used by the algorithm A n ^ t M- We have 



min(s,j^ j5 ) min(s,j^ 4 ) 



n = n m = n m ^ n i i+2 

i=i i=i 



log M 



1/6,. 



, . . j y ya^loga;- 1 



log oj 1 



, log M \ 5 ,. / log M 



log CO 1 



B+(log3)/03 5) 
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(?((! + log e-i^W*)) 



where the factor in the big O notation depends only on f3 and 5. This proves SPT+UEXP 
with 

^ log 3 

Since (3 can be arbitrarily close to one, and 6 can be arbitrarily close to a*, the exponent 
r* of SPT is at most 

^ log 3 

a* 

where for a* = oo we have 1°4^ — o. This completes the proof of Theorem [1] for the 
class A std . 



9 Relations to multivariate integration 

Multivariate integration 



INT S (/) 



f(x) da; 



[o,i] s 



for / from the Korobov space H{K s>aj b) was studied in [J]. It is easy to see that mul- 
tivariate approximation is not easier than multivariate integration, see e.g., [6]. More 
precisely, for any algorithm A ntS (f) = J2 k=1 (Xkf(xk) for multivariate approximation us- 
ing the nodes Xi, . . . , x n G [0, l) s and ot k G L 2 ([0, define (3 k '■— f^ Q ^ s oi k {x) dx and 
the algorithm 



k=i 



for multivariate integration. Then 
f(x)dx-A^(f) - 

[0,1]= 



\a k {x) f(x k ) dx 



< 



\ (/(*)-£ 

M&> V k=l 

I \f( x ) - y2 a k(x) f(x k )\ dx 

J tw s \ k=l J 



1/2 



= ||/-^n, S (/)|U 2 ([0,l] s )- 

This proves that for the worst-case error for integration we have 



sup 

feH(K s . a . b ) 
\<i 



„ n 

/ f(x)dx - y2/3 k f(x k 

Mi]' k=1 



<e L ^(H(K s ^ b ),A ntS ). 



e(H(K Statb ),A^ s ) 

f 

Since this holds for all algorithms A UiS we conclude that 

e int (n, s) := inf e(H(K s . a . b ), AjfJ < e app (n, s) := e 1 ^^ (n } s). (26) 
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Here e mt (n, s) and e app (n, s) are the nth minimal worst-case errors for multivariate in- 
tegration and approximation in H(K s>a ^), respectively. Furthermore for n = we have 
equality, 

e int (0,s) = e app (0,s) = 1. 
From these observations it follows that for e G (0, 1) and s G N we have 



n ie, s 



< n L2 - app ' Astd (e,s), (27) 



where n mt (e, s) is the information complexity for the integration problem. 

Obviously, for multivariate integration only the class A std makes sense. The inequal- 
ities fTSBl) and (|27p mean that all positive results for multivariate approximation and the 
class A std also hold for multivariate integration. In particular, the following facts hold: 

• Exponential convergence holds for multivariate integration for arbitrary a and b 
with the largest rate p mt (s) > 1/B(s). Although only uniform exponential conver- 
gence was considered in [I], the proof presented there allows to conclude that we 
have p mt {s) = l/B(s). 

• Uniform convergence holds for multivariate integration iff B < oo and the largest 
rate [p mt ]* = 1/B, as for multivariate approximation. This was shown in [I]. 

• Polynomial tractability and strong polynomial tractability for multivariate integra- 
tion were studied in [I], where it was shown that they are equivalent and hold iff 
B < oo and a/s are exponentially growing with j. These conditions are the same 
as for multivariate approximation. 

The exponent [r mt ]* of SPT for multivariate integration was estimated in [3], and 
we have [r mt ]* G [B,2B]. From Theorem [T] it follows that 

lyinti* < rappi* < g + 

L 1 ~ 1 ' ~ a* 
where [r app ]* is the exponent of SPT for the approximation problem. Hence we have 

log3 N 



[r mt ]* G 



B,B + min B 



Q 



which is an improvement of the result from [1] whenever a* > (log3)/i? which 
means that aj > exp(j'(a* — 5)) for large j. If a* = oo then 

j r intj* _ ^appj* _ g 

This is the case when dj > (1 + a) hj for large j and a > 0. 

Weak tractability for the integration problem was considered in [1] with a more 
demanding notion of WT. Suppose that we relax the notion of WT from [4], and 
use the notion of WT studied in this paper. That is, we say that the integration 
problem is weakly tractable if 

lim l0gr T M) = 0. (28) 

s+loge-i-s.oo S + log£: 1 
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We stress that the notion of WT as discussed in [I] implies (j2"8]l . but this does not 
hold the other way round. 

Using the definition (|28|) . we now show that we have the same condition linx, cij = oo 
for WT for the integration and approximation problems. Indeed, by Theorem HJ 
the condition linx, aj = oo implies WT for the approximation problem, which, by 
( |2"T|) . also implies WT for the integration problem. To show the converse, assume 
that the a/s are bounded, say aj < A < oo for all j G N. From ! 1 . Corollary 1] it 
follows that for all n < 2 s we have 

e int (n,s) > 2~ s ' 2 u 2 ~ 1 ^^ > 2~ s ' 2 uj m ' 2 = V s , 

where r\ := (co A /2) 1 ^ 2 G (0, 1). Hence, for e = r/ s /2 we have e mt (n, s) > e for all 
n < 2 s . This implies that n mt (e, s) > 2 s and 

logn int (e,s) slog2 log2 

; — T > --, 1 TT ~~ ^ — r > as S — y OO. 

s + log e 1 s + log 2 + s log ?7 1 1 + log r] 1 
Thus we do not have WT. 

This means that WT holds in the sense of f[28]) for the integration problem iff 
linij aj = oo, which is the same condition as for the approximation problem. 

Since for the integration problem we have UEXP iff B < oo, see jH Theorem 1], it 
follows that we have WT+UEXP iff B < oo and linx,- aj = oo. Again, this is the 
same condition as for the approximation problem. 



10 Analyticity of functions from H(K s a b 

In this section we show that the functions from the Korobov space H(K sa b) are analytic. 
Proposition 2 Functions f G H(K s a b) are analytic. 

Proof. Since H(K s a b ) C H(K s l l ) it suffices to show the assertion for / G H(K s l l ). 

Let ct = (ofi, ci2, . . • , a s ) G Nq with |a| = a.\ + • — h a s . For / G H(K s l ^), consider 
the operator D a of partial differentiation, 



cte" 1 dx 



Q2 
2 



Then 



7(h) (27ri)wn^ 



exp(27ri/i • £c) 



where, by convention, we take 0° = 1. 

Let 0J\ G (u, 1) and g = u/ui < 1. For any a G N consider ^(x) = x 2a q x for x > 0. 
Then g'(x) = if x = 2a/logg _1 and 



g" {2a/\ogq- 1 ) = \{^ 



2a 



a 



logg 



-i 



2a- 1 



logg < 0. 
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Hence, 



Since 



then 



g(x) < g (2a/ log q : ) 



2a 



e logg 



-i 



2a 



a 



2a 



.a 



al — - J <e 2a (a!) 2 
a\ 



g(x) < 



Hence, we have 



x 2a u x < 



\ogq- 1 



2a 



2<i 



fa!) 2 . 



log UJi — log CO 



(a!) V == C 2a (a!) V. 



Note that C depends only on u and u>i. 
Then tu h = implies 



\D a f{x)\ = 



< ll/l|tf(JC,l,l) 

< 11/11*^,1,1) 
< 
< 



1/2 



w i /2 (27ri)wn^ 

J'=l 



E(27r) 2H ni^| 2Q ^" jl 



exp(27ri/i • a;) 

1/2 



1/2 



||/||^,i,i)(27rC)HnK')(l + 



j'=i 



1 — o>i 



1/2 



s/2 



s/2 



where d = H/H^,,,,) (l + y^)' > and C 2 = 2nC > 0. 

Then for any £ = (d, . . . , ( s ) and any x = (x±, . . . , x s ) with ||a; — Clloo < C^ 1 we have 



v D " m f](x -cy 



aGNg i=l 
/ oo ^ 

< d ( X)(C 2 ||as - Clloo)*" 



v a=0 



= Cx 



1 — C 2 ||cc — CIU 



< oo. 



Hence / is analytic, as claimed. 



□ 
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